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Abstract 

We discuss the implementation of quasiparticle calculations for point de- 
fects on semiconductor surfaces and, as a specific example, present an ab 
initio study of the electronic structure of the As vacancy in the +1 charge 
state on the GaAs(llO) surface. The structural properties are calculated with 
the plane-wave pseudopotential method, and the quasiparticle energies are 
obtained from Hedin's GW approximation. Our calculations show that the 
la" vacancy state in the band gap is shifted from 0.06 to 0.65 eV above the 
valence-band maximum after the self-energy correction to the Kohn-Sham 
eigenvalues. The GW result is in close agreement with a recent surface pho- 
tovoltage imaging measurement. 

PACS: 71.15.-m; 71.45. Gm; 73.20.Hb 



1 Introduction 

Over the past 15 to 20 years intense research efforts have been made to understand 
the nature of native point defects near and at semiconductor surfaces. Such defects 
play an important role in surface electrical characteristics, e.g., Fermi-level pinning 
and charge-carrier recombination. They also influence the kinetics of adsorption, 
diffusion, and growth as well as surface chemical activity. The progress of the 
research efforts in understanding point defects on III-V semiconductor surfaces has 
recently been reviewed by Ebert jl]. 

At the atomistic level significant insight into the identification of these defects 
comes from scanning tunneling microscopy (STM). However, STM measures the 
local electronic density of states, which cannot always be directly interpreted as 
the atomic geometry. Therefore, accurate calculations have turned out to be very 
valuable for the interpretation of the experimental findings. In this paper we present 
an ab initio study of the arsenic vacancy Vas on the GaAs(llO) surface under p-type 
conditions. In STM images of the fiUed-state As sublattice, this vacancy gives rise 
to two distinct features: a dark hole of the lateral dimension of one As dangling 
bond and a long-range charge- related depression of the neighbouring As atoms |^ . 
The first feature is directly related to the missing atom, and the second one is due 
to a downward local band bending, from which it was concluded that in p-GaAs 
the arsenic vacancy is positively charged. STM images acquired under positive 
bias probe the empty pz orbitals of the Ga sublattice and show an enhancement 
of the contrast from the two Ga atoms nearest to the vacancy. Lengel et al. 
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interpreted their finding as an outward relaxation of the two Ga atoms neighbouring 
the vacancy. They also found support for this interpretation from a tight-binding 
calculation. Their calculation suggested a charge state of +2 for the vacancy. On 
the other hand, by comparing STM images of the arsenic vacancy with the images 
of Zn dopant atoms on p-type GaAs(llO), Chao et al. Q were able to determine 
the charge state of Vas to be +1, utilizing the fact that the dopant atom is in a 
charge state of —1 a priori. 

Recent density-functional theory (DFT) calculations by Zhang and Zunger Q 
and by Kim and Chelikowsky ||, ^ have confirmed the stability of the -1-1 charge 
state. Furthermore, their calculations have independently shown that the energy 
minimization of the defect geometry results in a downward relaxation of the neigh- 
bouring Ga atoms, and the calculated STM images show an enhancement of these 
two atoms in agreement with experiment. Both of these calculations were per- 
formed at a comparable level of sophistication but differed in the details. Zhang 
and Zunger [Q found an asymmetric relaxation of V^, whereas Kim and Che- 
likowsky 1^,^ found a symmetric relaxation, which seems to agree with experiment. 
However, in a recent combined experimental and theoretical study on a related 
system, the phosphorus vacancy on InP(llO), Ebert et al. suggested that the 
symmetric feature of the experimental STM image could also be explained by a 
thermal flip motion from two degenerate asymmetric configurations. 

For properties other than structural ones the experimental situation is less com- 
plete. There are so far no direct measurements of either defect ionization or charge- 
transfer levels for Vas on GaAs(llO). Indirect information about these quantities 
can be obtained, for instance, from measurements of the local band bending. Both 
the ionization and charge-transfer levels can be calculated theoretically, but so far 
this has not been done beyond DFT and using the local-density approximation 
(LDA). Furthermore, we note that the reported values deviate among the different 
studies. DFT-LDA also suffers from the underestimation of the fundamental band 
gap, which introduces a large uncertainty when interpreting the Kohn-Sham eigen- 
values of defect levels that lie inside the band gap as actual energy levels. As we 
will show in this paper, the band-gap problem for defect states can be circumvented 
by employing the GW approximation [|| for the electronic self-energy and calculat- 
ing the proper quasiparticle energies. The GW approximation has previously been 
applied to the perfect GaAs(llO) surface [p|, [To[| a s well as to bulk defects in some 
other materials, such as the Li vacancy in LiCl|l|] and the oxygen vacancy m zirco- 
nia , with good agreement between the theoretical predictions and the available 
experimental data. We note in passing that all the mentioned studies relied on 
plasmon-pole models, whereas our study takes the full frequency dependence of the 
screened Coulomb interaction into account. 

The paper is organized as follows. The method and computational details are 
described in Section]^ and results from the calculations in Section^. Finally, Section 
^ summarizes our conclusions. 



2 Method and Computational Details 

2.1 Atomic geometry 

The ground-state properties, such as the geometric structure of the defect, are 



calculated using DFT and the plane- wave pseudopotential method as imple- 

mented in the FHImd code |l^,|l6j. The exchange-correlation functional is treated 
in the local-density approximation as parametrized by Perdew and Zunger [ pTf . The 
single-particle orbitals are expanded in plane waves with a cutoff energy of 10 Ry 
in most of the calculations, but convergence tests were performed up to a cutoff 
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Fig. 1: A 4x2 surface unit cell with an As vacancy. The As atoms are shown in 
light grey and the Ga atoms in dark grey. The three Ga atoms next to the vacancy 
are indicated by arrows. Note the five-fold coordination of the Gan in the second 
layer. 



of 15 Ry. With the basis set used we find the bulk lattice constant to be 5.55 A 
(neglecting zero-point vibrations), which is in agreement with other calculations for 
III-V semiconductors P,p^, p^ and 1.8% smaller than the experimental value at 
room temperature of 5.65 A. ]20|. 

The As vacancy on GaAs(llO), whose geometry is shown in Fig. |l|, is described 
using the supercell method with periodic boundary conditions. The supercell slab 
consists of a 4x2 surface unit cell with six atomic layers, with one As atom missing 
in the top layer, and four vacuum layers. The dangling bonds at the bottom of 
the slab are passivated by pseudoatoms with noninteger nuclear charges of 0.75 and 
1.25 for As and Ga termination, respectively. We allow the atoms in the top three 
layers to fully relax while keeping the atoms in the three bottom layers at their 
theoretical bulk positions. In the case of the positively charged defect we apply a 
uniform compensating charge throughout the unit cell in order to maintain overall 
charge neutrality. The integration in reciprocal space was done using the special 
k-point (i, i, 0) |2^ in the irreducible part of the Brillouin zone. 

In order to test the influence of the supercell size we also performed a series of 
investigations using smaller 2x2 and 2x1 unit cells and found that the dispersion 
of the defect level, due to the spurious interaction of the defect with its periodic 
images, changed from 0.6 eV (2x2) to 1.5 eV (2x1), compared to 0.3 eV in the 
case of the 4x2 unit cell. For the 2x1 unit cell the defect-defect interaction is so 
strong that the two defect levels in the band gap exhibit a crossing, thus making 
any further calculations using such a small unit cell doubtful. On the other hand, 
the 2x2 unit cell gave results consistent with the ones obtained for the larger 4x2 
cell. In particular, the formation energy for the neutral defect turned out to be the 
same using either of the unit cells. 

Furthermore, we did a series of comparisons between six-layer slabs and four- 
layer slabs for the 2x2 unit cell and found no significant differences in the defect 
band structures (dispersion relations) for any of the possible charge states -|-1, 0, 
or —1. In the case of the four-layer slab only the top two layers were relaxed in the 
geometry optimization. 
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2.2 Quasiparticle energies 

Green-function theory is a mathematical framework for calculating the actual quasi- 
particle band structure, i.e., electron addition and removal energies, which charac- 
terize the excitation from an iV-particle system to an (iV±l)-particle system during 
a photoemission process. These energies, e^, are in principle given exactly by the 
quasiparticle equation 

(-^V^ + VUr) + Hj(r)) ^PZir) + J I](r, r'; CJ^^kCr') d'r' = eZ^Hir) , (1) 

where Vcxt is the external potential, Vn the Hartree potential, and ip'^^ the quasi- 
particle wave function. S denotes the self-energy operator, which contains all con- 
tributions from dynamic exchange and correlation. In general the self-energy is 
nonlocal, energy-dependent, and has a nonzero imaginary part, whose inverse is 
proportional to the lifetime of the excited state. For real systems S cannot be cal- 
culated exactly and must be approximated by a suitable functional expression. In 
the GW approximation E is given by ||^ 

E(r, r'; ir) = iG{r, r'; iT)W{r, r'; ir) , (2) 

where G is the one-particle Green function and W the dynamically screened Coulomb 
interaction. All operators are here written in real space and imaginary time. This 
representation has the advantage that the self-energy is a simple product, which is 
exploited in the GW space-time method that we have used for our calcula- 

tions. It avoids the multidimensional convolutions in reciprocal space and on the 
frequency axis that must be evaluated in conventional implementations. The inclu- 
sion of dynamic screening in the self-energy describes the correlation hole around 
individual electrons due to the Coulomb interaction. In this way the GW approx- 
imation goes beyond Hartree-Fock theory, which only describes exchange effects 
and ignores the correlation of electrons with different spin. Instead of simplified 
plasmon-pole models, which have been employed in almost all GW calculations for 
surfaces so far, we use the random-phase approximation with the full frequency 
dependence for the screened Coulomb interaction. At the end of the calculation 
the self-energy is Fourier transformed and analytically continued to the real energy 
axis. 

As usual, also in our calculations the Green function is obtained from the Kohn- 
Sham eigenfunctions, incorporating a large number of unoccupied states. The LDA 
wave functions form a reasonable starting point for the study of electronic band 
structures and are obtained from an equation closely resembling Eq. (P, with the 
local exchange-correlation potential, Vxc, in the place of the self-energy. The quasi- 
particle energies are therefore calculated from first-order perturbation theory ac- 
cording to 



Z = ^nk^ + (^^'^^^iSCerj - Vxcli^l:^^) , (3) 
where ej^^^ is the corresponding LDA eigenvalue associated with the wave function 

/LDA 
'Puk ■ 

Before investigating the system we first perform accurate convergence tests 
for the perfect surface, which can be represented with a smaller 1x1 surface unit 
cell. Fig. H shows the quasiparticle band structure for the GaAs(llO) surface. The 
dashed lines are obtained from DFT-LDA with a cutoff energy of 15 Ry, and the solid 
lines are the GW quasiparticle energies. The projected GW bulk band structure 
is indicated in grey. The ground-state density used for the DFT band-structure 
calculation was obtained using four special Monkhorst-Pack k-points |2^, which 
ensures the same BriouUin-zone sampling as a single k-point for the larger 2x2 unit 
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Fig. 2: Calculated surface band structure of relaxed GaAs(llO): dashed lines, DFT- 
LDA; solid lines, GW. The grey-shaded regions show the projected GW bulk bands. 



cell. The self-energy operator was evaluated at the four high-symmetry k-points 
(r,X,M,X'). The LDA surface bands are matched to the projected bulk bands by 
aligning the electrostatic potential in the central region of the surface slab with that 
of the bulk. The GW approximation has only a small effect on the dispersion of the 
occupied surface band and docs not change its relative distance to the valence-band 
edge at F. The GW and LDA occupied surface bands are hence aligned at this 
point. From Fig. || we then observe an almost uniform upward shift of the GW 
corrected unoccupied surface band by 0.8 eV, which is in close agreement with Zhu 
et al. Q and 0.1-0.2 eV less than found by Pulci et al. This opening of the gap 
makes for the typical improvement of GW quasiparticle band structures over DFT- 
LDA for the unoccupied states. We also note that the shift in the surface band is 
similar to the shift we find for the conduction bands in the bulk, 0.7 eV, illustrating 
the similarity of the wave-function character between the surface states and the 
bulk valence and conduction bands. The filled surface state at and in the vicinity of 
the F-point appears as a surface resonance, i.e., hybridization with extended bulk 
states is noticeable, and the peak becomes broad. The band structure at this point is 
shown for the sole purpose of illustrating the LDA and GW alignment. We repeated 
the calculations using different cutoff energies and found that the GW shift of 0.8 
eV remained for 13, 11, and 10 Ry cutoff. Also, for a four-layer slab we obtained 
a GW correction of 0.8 eV, and at 10 Ry cutoff we investigated the sensitivity 
of the GW corrected band structure to variations in the number of unoccupied 
bands. The self-energy is fully converged if 1049 unoccupied bands are included. 
The difference in the quasiparticle energies around the band gap is less than 0.05 
eV when the number of unoccupied bands is reduced to 153. With 379 unoccupied 
bands the agreement is better than 0.02 eV, which suffices for our purpose. For the 
GW calculation of the defect levels in the 2x2 surface unit cell we hence use 1500 
unoccupied bands, which correspond to 379 bands in the case of a 1 x 1 cell. 
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3 Results and Discussion 



3.1 Structural properties 

Let us first recall the structure of the perfect GaAs(llO) surface, which corresponds 
to the defect-free right-hand part of Fig. |^. This surface belongs to the Cih point 
group, which consists of a single mirror plane perpendicular to the [110] direc- 
tion. The relaxation of the surface preserves the point-group symmetry and consists 
chiefly of an inward movement of the Ga atom and an outward movement of the As 
atom. Hence the relaxation can be characterized mainly as a pure bond rotation, 
which results in a buckling of the surface fl^ . The Ga atoms rehybridize from an 
sp^ to an sp^ bonding situation to form a locally planar structure. The empty Pz- 
like orbital perpendicular to this plane forms the unoccupied surface band. The As 
atom maintains its sp"^ hybridization and binds to three Ga atoms in a pyramidal 
arrangement with a nonbonding electron pair in the fourth direction pointing away 
from the surface. The filled surface band is composed from the As lone electron 
pairs. 

We now turn the discussion to the As vacancy. The removal of the As atom 
leaves a dangling bond on each of the three Ga atoms surrounding the vacancy. 
As mentioned earlier, there is conflicting evidence for ^, |^ and against ||, ^ a 
breaking of the mirror symmetry due to the relaxation of the positively charged 
anion vacancy, which requires further careful studies. In this work we choose the 
second scenario and perform a relaxation of with the constraint that the Cih 
symmetry is preserved. The calculated equilibrium geometry shows a downward 
relaxation of the two Ga atoms in the first layer (Gai) and an upward movement of 
the second- layer Ga atom (Gan), which forms two new weak bonds with the Gai 
atoms. Thus the coordination of Gan is increased to five. In Fig. |l| the Gai and 
Gail atoms are indicated by arrows. We find the bond length between Gai and 
Gail to be 2.89 A for V^, which is slightly longer than the value of 2.78 A reported 
in Ref. ||]. When going from to V^, the Gai-Gan bond length decreases, 
indicating the bonding character of the la" defect level. 

3.2 Electronic properties 

The arsenic vacancy gives rise to three electronic states, la', la", and 2a', where a' 
denotes states that are even with respect to the mirror plane and a" denotes a state 
that is odd. The la' state is located several eV below the valence-band maximum 
and is thus always filled The la" state is located in the Kohn-Sham band gap, 
and depending on the location of the Fermi level, it is either empty or filled with 
one or two electrons. The la" state has a nodal plane coinciding with the mirror 
plane but has bonding character between Gai and Gan as mentioned above. The 
2a' state is found close to the conduction-band minimum (CBM) in the case of 
and clearly above the CBM for V^^ and V^s- 

Since the Kohn-Sham approach 
underestimates the band gap, it is plausible that the energies of the la" and 2a' 
vacancy levels are incorrectly given by DFT-LDA. We will focus on the la" level for 
y^. This state is mainly located at the Gai-Gan bonds and thus has a character 
different from any other in either the surface or the bulk, where there are no Ga-Ga 
bonds present. Therefore, one cannot a priori predict the quasiparticle correction 
from the shift of the normal surface or bulk bands; it has to be calculated separately. 
Our GW calculations for the shift of the defect level are done for a four-layer 2x2 
unit cell, using the convergence parameters that we found satisfactory for the clean 
surface. The defect-defect interaction gives rise to a rather large dispersion and 
requires a careful treatment in order to extract the most accurate value for the 
defect level. For this reason we fit the calculated DFT-LDA defect band to a simple 
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Fig. 3: Calculated band structure of V^^ on GaAs(llO): dashed lines, DFT-LDA; 
solid lines, GW . The horizontal lines mark the actual la" defect level as obtained 
at the special k-point. 

tight-binding (TB) model. Since the la" state is odd with respect to the mirror 
plane, one can, for the purpose of TB, regard it as an orbital with p symmetry that 
forms TT-type bonds in the [001] direction and exhibits a-bonding along the [110] 
direction. Within our supercell approach the defect sites form a rectangular lattice 
with lattice parameters and Cj,, where is the length of the supercell along 
[110] and ay the length along [001]. We found that we could fit the dispersion of 
the DFT-LDA (and also the GW) defect band by 

e(k) = ed + 2t/icr cos(fc2;a:c) -I- 2Vi^ cos(fcj,aj^) -f V2 cos(fc2,a2;) cos(fcj^aj^) 

+ 2V3CT cos(2fc^aa;) -I- 2V3^ cos{2kyay) , (4) 

where the fitting parameters ed and Vnrj have the meaning of the energy of a single 
defect and the hopping integrals, respectively. In order to properly reproduce the 
Kohn-Sham dispersion we included interactions up to third-nearest neighbours (Vso- 
and Vstt)- The correlation between the DFT-LDA results and the TB fit turned 
out to be 0.9999. At the special k-point kg = ^(a~^,a~^,0) the contribution 
from the first and second neighbours vanishes and e(ks) = ed + '^Vsa cos{2kxax) + 
2V3Tr cos{2kyay) . As the contribution from the third neighbours is very small, 0.06 
eV, we feel satisfied to directly take the DFT-LDA calculated value at kg as the 
defect level. Furthermore, the TB fit to the GW dispersion gives a similar estimate 
of Vsa and Vs^r as in DFT-LDA, so when we calculate the quasiparticle correction 
to the defect level the errors largely cancel. The results for are shown in Fig. 
H where the DFT-LDA results are marked with dashed lines and the GW results 
with solid lines. The GW and the DFT-LDA results are aligned at kg using the 
occupied surface band derived from the As electron lone pairs similarly as was done 
for the perfect surface. The As atoms contributing to this band are all situated in 
the defect- free row of surface As atoms along [110]. In Fig. ^ this state appears in 
the projected bulk bands (grey shaded) at —0.8 eV. From our calculation we find 
the quasiparticle correction to the la" state to be 0.59 eV. However, the four-layer 
slab is slightly too thin to allow an accurate alignment of the Kohn-Sham band 
structure to the projected bulk bands, so for this purpose we use a six-layer 4x2 
unit cell, for which we also expect the agreement between the actual defect level and 
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the value at kg to be even better than in the case of the 2x2 cell. At the DFT-LDA 
level the defect state is found to be 0.06 eV above the valence-band maximum and 
is indicated by a dashed horizontal line in Fig. ^. The GW correction adjusts this 
upward to 0.65 eV (solid horizontal line). The calculated values for the 2a' state 
are also shown in the figure. 

So far no direct photoemission measurement for this defect level has been re- 
ported in the literature. Hence we compare our results to a measurement by Aloni 
et al. Il^,^, who used atomically resolved surface photovoltage imaging with STM. 
They measured the tip-induced band bending by scanning along the [001] direction 
of the surface and found a pinning of the band bending at 0.53 eV, at the location 
of a single defect. When the tip bias was changed, the band bending remained the 
same at the defect site until the tip-induced band bending in the defect-free region 
exceeded 0.53 eV. At further negative tip-sample bias the signature of the defect 
disappeared, and the band bending was the same at the defect site as in the defect- 
free regions, from which one can conclude that at this point the tip-induced band 
bending had pushed the defect state clearly below the Fermi level, which changes 
the defect into the neutral charge state. From the knowledge of the bulk Fermi level 
the authors concluded that the defect level is located at 0.62±0.03 eV above the 
valence-band maximum for a flat-band situation. Our result of 0.65 eV turns out 
to be remarkably close to the experimental value, although it is not entirely clear 
to which extent the latter also includes energy contributions from atomic relax- 
ation processes. We estimate that these are small, however, because our DFT-LDA 
total-energy calculations for the and V^^ geometries, both with one electron 
occupying the defect level, indicate a maximum energy gain of only 0.16 eV for a 
complete relaxation. In any case, both the experiment and our GW calculation 
agree that the defect level is closer to mid-gap than to the valence-band maximum. 

4 Conclusions 

We have found the GW approximation to be a useful tool for investigating the 
electronic structure of point defects on semiconductor surfaces. The methodology 
and the possible difficulties are described in this paper. As an example we have 
studied the arsenic vacancy on GaAs(llO) in its positive charge state. While the 
LDA underestimates the fundamental band gap and places the unoccupied la" 
defect state just 0.06 eV above the valence-band maximum, the GW approximation 
opens the gap and shifts the defect level upward to 0.65 eV. This theoretical result 
is in good agreement with a recently reported experimental value of 0.62 eV that 
was deduced from atomically resolved surface photovoltage imaging with STM. 
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